Thermoelectric phenomena in disordered open quantum systems 
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Using a stochastic quantum approach, we study thermoelectric transport phenomena at low tem- 
peratures in disordered electrical systems connected to external baths. We discuss three different 
models of one-dimensional disordered electrons, namely the Anderson model of random on-site en- 
ergies, the random-dimer model and the random-hopping model - also relevant for random-spin 
models. We find that although the asymptotic behavior of transport in open systems is closely 
related to that in closed systems for these noninteracting models, the magnitude of thermoelectric 
transport strongly depends on the boundary conditions and the baths spectral properties. This 
shows the importance of employing theories of open quantum systems in the study of energy trans- 
port. 



Thermoelectric phenomena are attracting considerable 
attention both theoretically and experimentally, due to 
their fundamental unsolved aspects as well as their im- 
pact in energy conversion technology p}. Here, we study 
thermoelectricity in disordered systems and the role of 
localization Q, focusing mainly on energy transport by 
electrons. In real solids, electrical transport at finite tem- 
peratures is also controlled by inelastic electron-phonon 
scattering. Thus, we will only consider energy trans- 
port by electrons at low temperatures in the quantum 
regime. Quite often, despite their intrinsic open quan- 
tum system character, thermoelectric phenomena in dis- 
ordered electrical systems are described within the frame- 
work of closed quantum systems, namely systems that 
do not exchange particles and/or energy with the envi- 
ronment [1, HI]. In fact it has been shown that energy 
transport in open disordered harmonic chains is depen- 
dent not just on the properties of the system, but also 
on the baths connected to it [f| , as well as the boundary 
conditions @. Also there is a large class of driven dis- 
ordered lattice gas models of particles hopping on a lat- 
tice and interacting through hard-core exclusion where 
the nonequilibrium transport properties depend on the 
boundary conditions [7|, |8j. It is then natural to ask 
whether there are other fundamental properties of ther- 
moelectric phenomena that are not captured by closed- 
system theories and an open quantum system approach 
is necessary. 

To the best of our knowledge there are only a few 
microscopic studies (without assuming an explicit func- 
tional form of electrical conductivity) on thermoelectric 
transport properties in disordered closed systems d, [l(| • 
Recently, energy transport by electrons in disordered 
chains has been studied within an open quantum sys- 
tem approach [TP ], but this study does not address the 
thermodynamic limit of system size and thus cannot be 
directly compared with the results of closed system the- 
ories. We therefore study thermoelectric transport in a 
few models of open one-dimensional noninteracting dis- 
ordered systems. Our main result is that although the 
asymptotic nature of transport at low temperatures in 
open systems is closely related to that in closed systems, 
the magnitude of thermoelectric transport depends on 



the boundary conditions and the baths' properties. 

We employ a stochastic (Langevin) approach to inves- 
tigate steady-state charge and energy transport in dis- 
ordered noninteracting tight-binding lattices. This for- 
malism has been extensively applied to study thermal 
transport in classical disordered lattices (5 , [ijj as well as 
quantum electrical [l3| and phononic [l4| systems. Its 
main advantage is that one can explicitly consider the 
effect of baths and system-bath coupling. In this pa- 
per, we consider only non-interacting electrons. Interac- 
tions could be included by, e.g., working with a stochastic 
Schrodinger equation in the context of time-dependent 
current-density functional theory, as it has been recently 
suggested in Ref. [l5|. However, we leave this study for 
future investigations. 

We consider three different one-dimensional (ID) tight- 
binding models of disordered electrons, namely, the An- 
derson model of random on-site energy (diagonal dis- 
order) 0, the random dimer model (short-range corre- 
lated disorder) and the random hopping model (off- 
diagonal disorder) [T3|. All states in the ID independent 
random on-site energy Anderson model (AM) are expo- 
nentially localized for any strength of disorder [l8j]; thus 
there is no localization-delocalization transition in ID. 
Now, one can have a localization-delocalization transition 
in ID by introducing correlations (short or long range) in 
the random variables. The random dimer model (RDM) 
is the simplest example of that, where one or both of 
the two possible random on-site energies e a and tb are 
random in pairs. It has been shown that when both site 
energies appear in pairs there exist two real critical points 
with critical energies e a and £b if |e a — e&| < It. where t 
is the constant hopping strength between sites The 
random hopping model (RHM) is an example of a ID 
model where a delocalized state appears at the band cen- 
ter even without any correlation in the randomness [l7| . 
The last model has many common features with a wide 
class of random spin chains such as random XY spin 
chains. All of the above results for the different disor- 
dered models have been derived for closed systems. Here 
we are particularly interested to know how these results 
are affected by the coupling with baths and what conse- 
quences we should expect on thermoelectric phenomena. 
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The full system consists of a disordered wire of N sites 
and two infinite baths being connected to the wire at the 
two ends. The Hamiltonian is 
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-E/=i tmci+i + c i+i c i) + T,i=i e i c l c h 

"7i Ea=l ( C a 4+1 + c a+l4.) with i = L,R, and 

. -7i(cf f ci +4^), V« fl = - 7 k(cf f cat + 4cf). 
Here q, and denote operators on the wire, the 
left and the right baths, respectively. The Hamiltonian 
of the wire is denoted by Hw> that of the left (L) and 
the right (R) bath by TC R (with i = L,R); the tunneling 
Hamiltonian between the wire and the left (right) bath is 
Vy l/R (Vwr)- The tunneling from the disordered wire to 
the baths is controlled by the parameters 7^. The AM of 
diagonal disorder is defined by constant electron hopping, 
i.e., ti — t for 1 — 1,2, ..N — 1, and a distribution for the 
identically distributed independent random site energies 
Q. In the RDM t x = t for I = 1, 2, ..N - 1 but there are 
only two values of on-site energy ei = e a or eb which are 
assigned randomly in pairs with probabilities p and 1—p 
respectively. Finally, in the case of the RHM, e; = for 
I = 1,2, ..N, and the parameters ti are chosen from an 
independent random distribution. 

We assume that each bath is in equilibrium at a spec- 
ified temperature Ti and chemical potential fa (with 
i = L, R) before coupling it with the wire. The cou- 
pling of the baths with the wire introduces noise and dis- 
sipation in the wire. We apply the stochastic approach 
following Ref.[l3[ to derive steady-state charge and en- 
ergy current in the wire. Let us define j p and j u as the 
particle and the energy current density, respectively. 



3 



1 

2^ 



dui T 1n (uj) (Jl - }r) (2) 

r = ^ r duj hcu T ^^) - Jr) > (3) 



where T 1N = ^ 2 1 ' L 2 1 ' R 2 p l (lu) Pr (lu)\G 1n \ 2 /h\ G = 
Z~ l , Z = $ - t(uj) and 

®lrn = ( w - t) &l,m + 4^l,m-l + -T^l.m+l' (4) 

n n n 

Here <t is the Hamiltonian matrix of the disordered 
wire and G is the full Green's function of the wire coupled 
with the baths. The self-energy correction E, coming 
from the baths, is a N x N matrix whose only non-zero 
elements are En and T*nn- In the following we set fi = 
1 and discuss the case in which 7l = Jr = 7, "f' L = 
j' R = 7', which corresponds to symmetric coupling to the 
baths. In this case, £> R {o->) = EJ(cj) = E(u>) = 7 /2 gi i i(w) 
where gii(u>) is the single particle Green's function of the 
isolated bath at the first site a = 1: 
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Also in Eqs.JH), fi = l/{exp[{u - fa)/k B T 2 ] + 1} 
is the Fermi function of the ith bath and pi(u>) — 
— Im[g± t i(u))]/Tr is the local density of states at the first 
site (a = 1) on the ith bath. It can be shown that T\n(w) 
is the transmission coefficient of an electron from the left 
to the right bath through the disordered chain at en- 
ergy u>. Interestingly for 7=1 (ideal baths) and 7' = 1 
(ideal contacts) the above formulation merges with the 
Landauer theory of transport. 

We are interested in the system properties in the ther- 
modynamic limit. We thus extend a technique originally 
developed by Dhar 5] to determine steady-state ther- 
mal currents in a disordered harmonic chain connected 
to baths to the present case of disordered electrical open 
quantum systems. This technique is a generalization of 
the popular recursive Green's function method [10( to 
open systems. We are interested in finding the asymp- 
totic system size (N) dependence of the steady state (j p ) 
and (j u ), where (• • •) denotes average over disorder real- 
izations. Following [j| 0] we now separate out the wire 
and the bath contributions in Gin and write them ex- 
plicitly: 
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A n (lj) = D hN -E(w)(D 2 ,jv + D hN -i) +S 2 (w)I>a,JV-i 

where Ajv(cj) is the determinant of Z and Di m is the 
determinant of the sub-matrix of beginning with the 
Zth row and column, and ending with the mth row and 
column. We can numerically compute the elements D\^ m 
efficiently by taking the product of the 2x2 random 
matrices XJ: 
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After a little algebra we find 
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We compute (Tin(<^>)) numerically using Eqs.((5]l2i- The 
above formulation can be used to find steady-state cur- 
rents at finite bias, but in this paper we are inter- 
ested in the linear response regime at low tempera- 
tures, i.e., Afj, « fx, AT « T and k B T << fx with 
A^i = [i L - p R , AT = T L - Tr, /i = (p L + fX R )/2 and 
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FIG. 1: (Color online) Plot of (Tin (-0.5)) versus N for three 
different sets of 7 and 7' in the AM with t = 1. at each 
site I is drawn from a Gaussian distribution of mean zero and 
variance a 2 — 0.0025. The full curve is asymptotic scaling of 
(71jv(— 0.5)) for all three sets of parameters. 



FIG. 2: (Color online) Plot of (T 1N (0)) versus N in the pinned 
and unpinned RHM for different 7, 7' with t = 0.8 and At = 
0.2. The two straight lines correspond to the asymptotic of 
(7ijv(0)). The sample to sample fluctuations in Tljv(0) are 
quite large (on the order of (Tljv(0))) and do not decay with 
increasing realizations or system size. 



T = (Tl + Tr)/2. In the linear response regime we write 
the heat current as j h = j u — ixj v . We then find 

Ztt o 

Ztt o o 

Results— We thus see from the above results that par- 
ticle and thermal conductances scale similarly even in the 
disordered open quantum systems and the Weidemann- 
Franz relation is valid in the absence of interactions. We 
also note that the asymptotic size dependence of the elec- 
trical and thermal conductances depend on the asymp- 
totic behavior of (Tin(h))- Hereafter we will then focus 
on this quantity. We find that for independent random 
distribution of on-site disorder in the ID AM, both parti- 
cle current and heat current become exponentially small 
in the asymptotic limit of system size in the closed as 
well as the open systems. Thus the behavior of (7ijv(m)) 
in the asymptotic regime is dominated by localization 
physics and there is clearly no diffusive behavior satisfy- 
ing either Ohm's or Fourier's law. The asymptotic char- 
acter of the eigenstates of disordered systems is quanti- 
fied through the localization length £(w) which is defined 
by 

IhnN^N-hn^N^)) = -2/e(w). (8) 

For a disordered tight-binding chain of lattice constant 
a, with a random on-site potential e/ at each site I drawn 
from a Gaussian distribution of mean zero and variance 
<7 2 , the localization length of the closed chain is given by 
£(E F ) = 2(a/a 2 ){At 2 - E 2 F ), where E F is the chemical 
potential at zero temperature or the Fermi energy [2l| . 
Thus i(E F ) = 3000 for a 2 = 0.0025 t 2 and E F = -0.5* 
with a = 1. We compute numerically the localization 



length in the disordered closed and open chains for the 
same above parameters using the definition of Eq.©. 
We find in numerics with t = 1, a = 1 (see FigJT]) that 
£,(E F ) (~ 6000) is the same for both the closed and the 
open chains (7' = 0.7, 7 = 1.2). We also find that the 
magnitude of (Tin(jj,)) in the open chains falls rapidly 
from that of the closed chains for any changes of 7 and 7' 
from the unity. This is expected and can be understood 
physically. Any value of 7 and 7' different from unity 
introduces extra scattering in the chain and thus reduces 
the magnitude of (Tin((x)). 

It has been shown in Ref.[lj| that all states of the 
closed RDM are localized except at the two energies e a 
and Cb which are real critical points with infinite £ if 
|e« — £6 1 < 2i. In numerics with the closed systems we 
instead find that though (Tin(lo)) at uj = e a or re- 
mains constant with increasing system size, it decays 
algebraically ((7ijv(o;)) <~ A^ -2 ) with system size for 
\e a — e&| ~ 2t. Thus it is hard to conclude convincingly 
whether the energies (i.e., u> = e a or tb) are true criti- 
cal points for \e a — e&| ~ 2t. Interestingly, the sample to 
sample fluctuations in (Ti^(u>)) for the latter case decay 
with increasing system sizes; while away from these ener- 
gies (or in the AM) the sample to sample fluctuations in 
(Tin(oj)) do not fall with increasing length. In the open 
systems (with 7, 7' ^ 1) we again find similar asymp- 
totic behavior as the closed systems in numerics. Thus 
in open systems the transport is ballistic at uj = e a or 
for \e a — eft I < 2t and it shows power law dependence on 
length at uj — e a or e& for \e a — e&| ~ 2t. 

It has been argued from simple considerations that for 
the closed RHM, a state at the band center, uj = 0, is 
extended and all other states are exponentially localized 
[llj ■ Here, again we find in numerics that (7ijv(0)) does 
not remain constant as a function of N in the closed 
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and open RHM, but it decays with increasing N. The 
asymptotic length dependence of (7ijv(0)) is given by 
jy-o.5 ^ gee pig(2]). We use a uniform distribution be- 
tween t and t + At for The RHM model is equiv- 
alent to a disordered linear chain of harmonic oscilla- 
tors for closed systems (T3, [HJ. Here, we wish to com- 
pare the open system results of the RHM with that of 
the random spring harmonic chains. Recently the au- 
thors of [HI have investigated heat conduction in ran- 
dom spring quantum harmonic chains for shorter lengths; 
but they were not able to conclude about the length de- 
pendence of the disorder averaged steady-state thermal 
current ((J)) in the spring model. It has been already 
argued in Ref. [f| that the asymptotic length dependence 
of the classical and quantum thermal currents are simi- 
lar in disordered harmonic chains. Also the asymptotic 
length dependence of (J) in the random spring harmonic 
chains, for two different models of baths (Rubin's baths: 

E(w) = k{\ - muj 2 /2k - lujim/k) 1 / 2 [l - mw 2 /(4/c)] 1/2 } 
with k the spring constant and m the mass of the lattice 
site; and white-noise baths: = —iju) are similar to 

that of the random mass harmonic chains: (J) ~ jV~°' D 
for Rubin's baths and (J) ~ iV -1 " 5 for white noise baths. 
Interestingly, we see that the asymptotic length depen- 
dence of thermal currents in the Rubin's bath case is sim- 
ilar to that of the open RHM. We further put two pinning 
potentials (ei = ejy = 1) at the two ends of the RHM to 
examine the effect of different boundary conditions. We 
find that the scaling of (7ijv(0)) remains the same as 
the unpinned case (see Fig(2]) . Two external quadratic 
pinning potentials at the two ends of the disordered har- 
monic chain with Rubin's baths change the asymptotic 
length dependence of (J) to iV -1 " 5 , thus showing a dif- 
ference in energy transport by the RHM and the ran- 
dom spring harmonic chain. This can be understood as 
follows. While energy transport in tight-binding chains 



mostly occurs by electrons at the chemical potential, the 
full band of conducting modes in harmonic chains carries 
energy. The external pinning does not affect the band 
center in the RHM; but breaks the translational invari- 
ance in the harmonic chains and pinches off the band 
of conducting modes from the zero frequency side, thus 
reducing (J). 

In conclusion, we have shown that the asymptotic na- 
ture of thermoelectric transport in the noninteracting dis- 
ordered open systems is quite similar to that in the closed 
systems. However, the magnitude of thermal and electri- 
cal conductances is smaller in the open systems compared 
to that in the closed systems. In earlier studies the ef- 
fect of coupling with baths has been included through 
a phenomenological lifetime due to inelastic scattering 
from the baths. This is done by energy continuation into 
the complex plane. Here, we have explicitly included the 
baths in our microscopic analysis. Since the technique 
used in this paper can be extended to higher dimensions 
following Ref. [24j , it would be interesting to analyze our 
results in 2D and 3D. It is known that the Anderson 
localization-delocalization transition in the random 3D 
AM will be smoothed in the presence of inelastic scatter- 
ing (due to baths), but it will be interesting to check how 
the corresponding transport properties will be affected by 
explicit coupling with the baths near the transition. Ex- 
periments in disordered systems are mostly carried out 
in open configurations. In fact, many real disordered 
systems such as doped polyaniline, random semiconduc- 
tor superlattices [25[ and random antiferromagnetic spin 
chains [26[ are considered to have similarity with the 
RDM and the RHM. Therefore, we expect our results to 
be useful in understanding experiments in these systems. 
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